Engineering Giant Nonlinearities in Quantum Nano-Systems 
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We describe a method to engineer giant nonlinearities in, and probes to measure nonlinear ob- 
servables of, mesoscopic quantum resonators. This involves tailoring the Hamiltonian of a simple 
auxiliary system perturbatively coupled to the resonator, and has the potential to engineer a wide 
range of nonlinearities to high accuracy. We give a number of explicit examples, including a readily 
realizable two-qubit auxiliary system that creates an potential and a x*"'' (Kerr) nonlinearity, 
valid to fifth-order in the perturbative coupling. 
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In the last few years there has been rapid progress in 
the development of nano-mechanical resonators and su- 
perconducting microwave oscillators [H, 0, S Both 
kinds of oscillators can be combined with other electri- 
cal elements to form mesoscopic devices, and because of 
their high frequencies, both have the potential to realize 
quantum behavior. If either kind could be built with non- 
linearities that were strong in the quantum regime (so- 
called giant nonlinearities [5|), this would open up many 
potential applications. These include quantum comput- 
ing simulating many-body systems 0], and studying 
the quantum-to-classical transition However, directly 
constructing resonators with such nonlinearities is not 
possible, at least with current technology 0. 

Measuring nonlinear observables of resonators is also 
important; it reveals signatures of quantum dynam- 
ics [13, 111, and we expect it to have future applications 
in feedback control and adaptive measurement [l^, [13 ■ 
To do so one must find mesoscopic devices that have in- 
teractions with resonators that are proportional to these 
observables. Physically building such interactions is a 
very challenging task. The only method suggested to 
date that could generate a range of nonlinearities or non- 
linear interactions requires active control of an auxiliary 
system on fast time-scales [3, 15 1. Here we present a 
method to do this that does not suffer this limitation. 

We recall first that it is well-known in the field of 
quantum optics that a low-dimensional system can gen- 
erate an effective nonlinearity in a larger system. For 
example, four-level atoms with the right level structures, 
when detuned appropriately from a cavity mode, will cre- 
ate x^^^ and other nonlinearities for the mode [^, [l^ . 
This is promising, because there are a number of low- 
dimensional mesoscopic systems that can be interfaced 
with resonators. However, to make this practical one 
must find systems with the right internal structure to 
generate the desired perturbation, and that can be read- 
ily realized. The simple and physically intuitive examples 
explored in quantum optics are natural for multi-level 



atoms, but not necessarily for mesoscopic systems. In 
addition, we wish to engineer nonlinear probes as well as 
nonlinarities. To achieve this goal we introduce a sys- 
tematic method for finding auxiliary systems to generate 
a given dynamics, including calculating this dynamics to 
high-order in the perturbation. This allows us to search 
over many auxiliary systems to find the simplest systems 
that provide a given nonlinearity to a given accuracy. 
The result is a method that can engineer both nonlinear- 
ities and nonlinear probes, and suggests that both the 
accuracy and the range of these nonlinearities can be 
increased by increasing the number of qubits in the aux- 
iliary system. (We note that this work was not, in fact, 
inspired by nonlinear optics, but by the recent work on 
"quantum gadgets" , perturbative auxiliary systems, in- 
troduced to connect other systems together, for enabling 
analyses of complexity in quantum computing 17|.) 

The basic idea is as follows. We couple a mesoscopic 
resonator, with annihilation operator a and dimension- 
less position X = [a-\- a^)/^/2, to a low-dimesional auxil- 
iary system via the linear interaction i/jnt = inxV . Here 
V is an operator of the auxiliary system, and /i gives 
the interaction strength. Interactions proportional to 
X are straightforward to engineer between mesoscopic 
resonators and charged systems such as Cooper-pair 
boxes [l3l), polar molecules [I^, or quantum dots [20| . 
How we design the system so that V is simple will be 
described below. Wc now choose the auxiliary system so 
that the separation between its adjacent eigenstates. A, 
is significantly larger than the maximum eigenvalue of 
the operator ^lxV that will be explored during the evolu- 
tion of the resonator. The joint Hamiltonian of the two 
systems may be written as 



H = Ho+ fixV + H,^ 



(1) 



where Hq is the Hamiltonian of the auxiliary system and 
fixV is a perturbation to this Hamiltonian. Since x com- 
mutes with all operators of the auxiliary, we can use stan- 
dard time-independent perturbation theory to diagonal- 
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ize i?aux = Hq + fixV. This gives us the eigenvalues 
of i?aux, E„, as a power series in A = fix. That is, 
_E„ = A ^^^g i!ji™^(/z/A)™a;™ for some real numbers 
Ek'^l If we now place the auxiliary system in one of the 
eigenvalues of -ffaux, say Ei, then the Hamiltonian for the 
resonator becomes ffong = A ^^^^^ iji™^£'"x"' + i?rcs, 
where e = /i/A. By designing Hq and F so as to elimi- 
nate specific terms in this sum, we will be able to generate 
nonlinearities of the form x™ for the resonator. 

For the auxiliary system to act as a probe to measure a 
nonlinear observable A, we must create an effective inter- 
action Hint = fJ,AV' where V' is an operator of the aux- 
iliary. Information about the resonator operator A can 
now be obtained by measuring an observable of the aux- 
iliary that does not commute with V' (this is merely the 
standard von Neumann measurement prescription [2l|). 
We can engineer a nonlinear interaction by tailoring the 
expansions of two eigenstates of the auxiliary, and placing 
the auxiliary in the space spanned by these eigenstates. 
In this subspace, the resulting Hamiltonian has the same 
form as i?eng above, but the numbers E'i'"' become op- 
erators in this subspace. By eliminating specific terms in 
this sum, we can obtain an interaction proportional to a 
given power of the observable A. 

Our task is not merely to obtain an ifo and V for the 
auxiliary that gives the desired perturbation expansion, 
but to find those that are simple to implement. To do 
this we note that a physical interaction fixV is usually 
straightforward to construct (i.e. natural) when V is di- 
agonal in the charge basis of the auxiliary system. We 
start be specifying a diagonal auxiliary Hamiltonian iJo, 
and perform the perturbation calculation to determine 
the required interaction operator V . We then find the 
unitary, [/, that diagonalizes V , and this gives us our 
diagonal V {V = UVU'^). The Hamiltonian for the aux- 
iliary system is then given by = UHqU^ ■ This is not 
diagonal in the charge basis, and is the Hamiltonian that 
we must engineer for the auxiliary system. For a given 
problem we will search for those auxiliary Hamiltonians 
that are the simplest to physically construct. 

In what follows we will consider using auxiliary systems 
of dimension two, three and four to generate nonlineari- 
ties. Note that to generate a nonlinearity z", we should 
eliminate from the perturbation expansion all terms of 
order less than n, since these will otherwise dominate. 
It thus makes sense to choose the diagonal elements of 
V to be zero to eliminate all first-order terms. It is also 
desirable to eliminate as many terms as possible whose 
order is higher that n. While these terms decrease with 
increasing order, to give an accurate rendering of cc" we 
will eliminate the term of order n + 1, and minimize the 
higher-order terms. To proceed we need the expressions 
for the coefficients of the various orders, £'i™\ of the 
perturbation expansions of the eigenvalues. We now give 
these up to fourth-order, where each expression for 



has been simplified by the condition that En = for 
1 < j < m. With this simplification the coefficients are 

Ei? = E I^^IVA™, (2) 

(3) 



= V;mKniV^«/(A„„Az„), 
p(4) — \ ^ VnmVmlVljVjn | Krin P I ^hi P /.n 

" ~ ^ A, A- A2 ' ^ ' 

, . , ^mn'—^L7i^in , , ^mn^]n 

m^L,j^n m.l^n 

where the Vnm are the matrix elements of V . Higher- 
order terms are significantly more complex, but are easily 
calculated numerically. 

To begin we consider a single auxiliary qubit. For a 
single qubit we are essentially restricted to H^ = {/S./2)(Jx 
and V = Uz- The resulting expansion has only terms of 
even order. To fifth order the effective Hamiltonian is 



i/eff = Act; [(1/2) - 2e^x^ + iie^x'^] + H,^ 



(5) 



where cr^ = (1 - 2e^ + &e^)a^ + (2e - Ae^)az We 
cannot engineer the nonlinear potential a;^ by itself with 
a single qubit, since we do not have the freedom to elimi- 
nate the second-order term without removing the fourth- 
order term. However, we see that a single qubit generates 
the interaction cr'^x^ accurate to third-order in e. This 
allows us to measure x"^ by using a sing le-electron tran- 
sistor to measure CTz for the qubit [22|. It was shown 
recently that this nonlinear measurement is useful, as it 
will generate mesoscopic-superpositon states of the res- 
onator directly from a thermal state [isj . 

In the above case, by choosing the resonator frequency, 
to be much larger than the strength of the interaction, 
/i = Ae, and making the rotating- wave approximation, 
we obtain the familar energy-coupling a'^a^a [23| . If the 
resonator frequency is much larger than the interaction 
strength, and we do not wish the coupling to degen- 
erate to the energy coupling, then we can modulate the 
interaction strength /i at a frequency, v, close to to. This 
neat trick effectively brings the resonator frequency down 
to Lo ~ v from the point of view of the coupling 11 1 . 

The above analysis shows that to generate x^, the 
auxiliary system must have at least three levels. We 
now examine what can be achieved with a three-level 
(qutrit) auxiliary. We take the energy levels of Hq to be 
evenly separated by A, so that En^ = nA, n = 0,1,2, 
and examine the perturbation expansion for the central 
level. El . From Eq. ^ we find that we can eliminate the 
second-order term by choosing V so that |Vi2p = l^oip. 
Conveniently, for a qutrit, this also removes the fourth- 
order term. Thus, placing the qutrit in the eigenstate 
corresponding to Ei, this single condition generates the 
resonator Hamiltonian 



i/e„g= (2ARe[FoiVl2V^2o]e' 



(6) 
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accurate to fourth order in e. Note that the condition 
I Vi2 P = I Vbi P gives us considerable scope in choosing the 
physical interaction operator V. In particular, we have 
complete freedom in choosing the diagonal elements of 
V. On the contrary, all the elements of Hq, the physical 
Hamiltonian for the qutrit, are non-zero. This means 
that one must couple all three charge states of the qutrit. 
If the qutrit is a polar molecule or quantum dot, then this 
could be achieved using three lasers tuned respectively to 
the three transitions. 

An auxiliary qutrit docs not have enough degrees of 
freedom to generate the nonlinear potential (while 
removing the second and third-order terms), or indeed 
to realize a probe for measuring x^. For these tasks 
one must use an auxiliary with at least four-levels (two 
qubits). We now turn to this case, and consider first en- 
gineering an potential. Once again we take the energy 
levels of Hq to be = riA, and examine the expansion 
of El. Using Eqs.([2l) ©i '^^ fi^^d that the condition 
required to eliminate the second-order term is |Vi3p = 
2(|Vbi ^ — |Vi2p), and to eliminate the third-order term 

is Re 21/10^02^^21 + VwVosVsi - ^121^23^^31 = . These 
two conditions also remove the fifth-order term. 

To find practical auxiliary systems, we perform a nu- 
merical search over all V that satisfy these conditions, 
minimizing the distance between the resulting auxiliary 
Hamiltonian, Hq, and those that can be readily engi- 
neered. So long as the two conditions are enforced, the 
optimization is not difficult and gives many solutions. A 
nice example is the two-qubit Hamiltonian 



Ho = A aa^y^af'' +hof^ -ca^^^'af^ 



(7) 



where a = 0.914, h = 0.405 and c = 0.5. The interac- 
tion terms in this Hamiltonian can be engineered between 
superconducting qubits using, for example, methods de- 
vised by Mooij tt al. . The required interaction with 
the target resonator is 



(1) 



(2) 



(8) 



where / — —1.823 and g = —1.382. This auxiliary sys- 
tem generates an expansion for Ei in which all the terms 
of odd-order are eliminated (at least up to seventh order) 
and the coefficients of the even terms are, E^^^ ~ 2x 10^"*, 



E^^^ = -1 and E^^^ = -3.99. It thus generates an 
potential accurate to fifth order in e. 

The above auxiliary system also allows us to gener- 
ate the Kerr (or x^'^^) nonlinearity, whose Hamiltonian 
is Hkctt oc (a^a)^. To do this we use exactly the same 
adiabatic elimination process that turns the interaction 
axx'^ into axa^a: we choose the resonator frequency to 
be much larger than the strength of the x'^ potential and 
make the rotating wave approximation. This eliminates 
all terms in x'^ that are not powers of a^a, transforming 
x^ into (3/2)(ata)2, and thus x^^^- 



p(6) 



Now consider engineering a probe for measuring (or 
equivalently (a^a)^). To do this we must tailor the ex- 
pansions of two eigenvalues. If we choose Ei and E2, then 
this merely involves using the conditions above for tailor- 
ing El , and adding to them the equivalent conditions for 
E2 (these can be obtained from the former by symme- 
try). Performing a numerical search to find Hq and V, 
our results suggest that there is only one solution, being 



Ho 



i'' + (1/2)4^' 



(9) 



.(1) 



(2)1 



with the interaction operator V ~ l^[f<^z 
where / = 1.682 and g = 1.189. This Hamiltonian 
is even simpler than the one presented in Eq.([71), as 
no coupling is required between the two qubits. This 
Hamiltonian can, of course, also be used to engineer x^ 
and x*-^-* nonlinearities by placing the auxiliary in the 
eigenstate with eigenvalue Ei (or E2). This configura- 
tion eliminates all the odd terms in the expansion for 
El and E2 up to seventh order, and generates the even 



terms E', 



(2) 



-E- 



(2) 



8 X 10" 



E. 



(4) 



-E- 



(4) 
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and i?2^' = —E[^^ = —4.24. Denoting the eigenstates of 
the auxiliary with eigenvalues Ei and E2 as |1) and |2), 
this auxiliary system generates the effective interaction 
Heft = As'^Yx'^, where Y = |1)(1| - |2)(2|. 

To show that our method can be realized with cur- 
rent technology, we consider implementing an x^ nonlin- 
earity using the configuration given in Eq.([5|). This is 
achieved by coupling the resonator to two Cooper-pair 
boxes (CPBs) via the usual capacitative interaction [25|. 
Operating each CPB at its degeneracy point gives ex- 
actly the right Hamiltonian; the Hamiltonian for the 



f'' CPB is Hj = 
Josephson energy |1 
(1,1/2) X lO^s-i 




where hojj is the (tunable) 
Reasonable values are (a;j',a;2) 
Scahng the resonator's position 



so that X — {a + a^)/^/2, current interaction strengths 
for superconducting resonators are fj, ~ 10® [3|, a nd this 
is also realistic for nanomcchanical resonators [25[ . With 
these parameters, and placing the qubits in the eigenstate 
of E2 (achieved using simple single-qubit operations) the 
engineered Hamiltonian is 



Heng = hK [x^ - 0.042 x^ + 0{e^)] + iJ,c 



(10) 



with K = 10^ s~^. This is a giant nonlinearity as 
advertised, albeit with a (relatively) small additional x^ 
term. With a resonator frequency of 100 MHz, this 
term is to excellent approximation a x^^' nonlinearity, 
given by the Hamiltonian Hcng = —{3/2)hK{a'^a)'^. As a 
benchmark for the strength of this nonlinearity, one can 
use the time it takes to create a mesoscopic-superposition 
of two coherent states, starting from a coherent state. 
This time is conveniently independent of the amplitude 
of the initial coherent state. The nonlinearity above per- 
forms this task in t = tt /(3k) = lO^s [3. (To compare, 
resonator damping times are ~ 10 ms for mechanical [l[, 
and 10/is for superconducting [sj.) 
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We now consider the robustness of the method to er- 
rors in the auxiUary Hamiltonian. To do this we add 
an independent Gaussian error with standard deviation 
O-OlW]' (errors of ~ 1%) to the three degrees of freedom 
of each CPB Hamiltonian, and calculate the resulting 
perturbation expansion for i?2- We then calculate the 
standard deviation of the induced change in each term 
in the expansion of E2 by averaging over many samples. 
Denoting the standard deviation of the change in the co- 
efficient £'2™'' as we find that these are very small: 
(S(i),...,I](7)) = 10^3 X (0.2,0.3,0.8,1-3,3.3,5.8,15.3). 
The method is thus quite robust to errors. This calcu- 
lation also provides an estimate of the errors induced by 
dcphasing, since this noise can be modeled as fluctuations 
in the Hamiltonian. Errors from dcphasing should be sig- 
nificantly smaller than those calculated above, since the 
best dcphasing rates for CPB's are currently 10^ s~^ ~ 
10~^uji [2^. We can also estimate the effects of noise due 
to decay of the upper levels of the CPBs in the same way. 
The broadening of the energy levels for each CPB is given 
by the non-Hcrmitian Hamiltonian H ~ — i?i(7/2)a+cr_, 
where 7 « lO'^uif is the decay rate. Including this 
Hamiltonian and re-calculating the perturbation expan- 
sion reveals that the changes induced in each coefhcicnt 
are no greater than 10~^. 

Returning to Eg. pU]) . we can reduce the relative size 
of the sixth-order correction by increasing A with respect 
to /i, which is quite practical. However, this will also re- 
duce the strength of the nonlinearity. If we want to 
both increase the cc^ term and reduce the term, then 
we need to reduce the coefficient of the sixth-order term 
in the perturbation expansion. Our numerical search in- 
dicates that there is a limit to how small one can make 
the sixth-order term with a 4-level auxiliary. We expect 
that using an auxiliary with 3 qubits (8 levels) will allow 
elimination of the sixth-order term completely. 

As useful as the results we have presented here are for 
engineering mesoscopic systems, this technique appears 
to promise even greater potential with the use of larger 
auxiliary systems (e.g. three qubits). Realizing this po- 
tential will require exploring the complex landscape that 
maps Hamiltonians to perturbation expansions, and this 
is an interesting question for future work. 
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